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We have studied the low-lying excitations of a chain of coupled circuit-QED systems in the ultrastrong cou- 
pling regime, and report several intriguing properties of its two nearly degenerate ground states. The ground 
states are Schrodinger cat states at a truly large scale, involving maximal entanglement between the resonators 
and the qubits, and are mathematically equivalent to Majorana bound states. With a suitable design of physical 
qubits, they are protected against local fluctuations and constitute a non-local qubit. Further, they can be probed 
and manipulated coherently by attaching an empty resonator to one end of the circuit-QED chain. 



Confronted with formidable difficulties in solving strongly 
interacting many-body systems, it has been desired to find 
good quantum simulators. It may seem natural to simulate 
a many-body system with another tunable system of massive 
particles such as ultracold atomic gases |1 1. In fact, any con- 
trollable quantum system, notably quantum computer if ever 
practical, can simulate efficiently many -body systems [2|. In- 
deed it has been recognized that photons confined in coupled- 
cavities simulate closely the quantum behaviors of strongly- 
coiTelated many -body systems |[3] |4l. Subsequent studies 
have revealed that Bose-Hubbard model f5 1, interacting spin 
models |4i i6J, and other exotic quantum phases [7| can be 
simulated efficiently using the coupled-cavities. Further, re- 
cent advances in solid-state devices such as circuit-QED sys- 
tems [8 9 1 and micro-cavities ifToHTTl and ongoing efforts to 
fabricate large-scale cavity arrays lfT2ll make the array of cou- 
pled cavities a promising candidate for an efficient quantum 
simulator 

Meanwhile, the ultrastrong coupling regime of the cavity- 
QED system, where the light-matter coupling energy is com- 
parable to or even higher than the energy of the cavity 
field, has been envisioned |fT3]| and experimentally demon- 
strated lfT4l . The ultrastrong coupling brings about fundamen- 
tally different physics deeply connected to the high degree 
of entanglement between the "matter" and the photon 1.151 - 
[T9l . However, the effect of ultrastrong coupling on the low- 
energy excitations of an array of coupled cavity-QED systems 
remains unclear, and is our main concern in this work. 

In this paper, we investigate the low-lying excitations of a 
one-dimensional (ID) array of circuit-QED systems (cQEDs), 
with each cQED being in the ultrastrong coupling regime; see 
Fig.[T] It turns out that the array permits two nearly degenerate 
ground states separated by a finite energy gap from the contin- 
uum of higher-energy states. We find several intriguing prop- 
erties of the two ground states: (i) They are Schrodinger cat 
states at a truly large scale, and involve maximal entanglement 
between the resonators and the qubits. (ii) With a suitable de- 
sign of physical qubits, the two ground states are protected 
against local fluctuations and constitute a non-local qubit ll20l . 
(iii) They are mathematically equivalent to the long-searched 
Majorana bound states ||2T|. (iv) They can be probed and ma- 
nipulated coherently by attaching an empty resonator to one 
end of the circuit-QED chain. Such configuration turns the 



total system (the circuit-QED chain plus the empty resonator) 
into another effective circuit-QED system. There are many 
promising types of superconducting qubits, among which we 
focus on Fluxonium |22|. As we illustrate below, its strong in- 
ductive coupling with the superconducting resonator 1 17 1 and 
its anisotropic noise characteristics [18J are well suited for our 
purpose. 

System: a circuit-QED chain — We consider a ID ar- 
ray of cQEDs; see Fig. [T] Each cQED consists of the "res- 
onator", a superconducting microwave transmission line, and 
the "qubit", a superconducting quantum bit (two-level sys- 
tem) |l8]|23l, and is theoretically described by the Rabi Hamil- 
tonian 
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where and a\ are the field operators of the resonator with 
frequency wo, erf and af Pauli operators of the qubit with 
energy splitting fi, and A the resonator-qubit coupling energy 
in the ith cQED. The resonators of neighboring cQEDs are 
coupled capacitively to each other, and photons hop from one 
resonator to nearby ones. The Hamiltonian of the whole chain 
is thus given by 
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FIG. 1. Schematic of ID circuit-QED arrays. The red dots indi- 
cate qubits placed inside of superconducting resonator. The TVth res- 
onator is coupled to the detection resonator. The circuit-QED array 
realizes the transverse field Ising model (TFIM), and the detection 
resonator can measure and control the degenerate ground state of the 
TFIM. 
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where J is the photon hopping amplitude and N is the number 
of cQEDs in the chain. 

Before discussing the energy levels and associated wave- 
functions of the whole chain, we first briefly review the prop- 
erties of the low-lying states of a single cQED in the ultra- 
strong coupling regime (A > luq). The strong coupling dis- 
ables the standard rotating wave approximation, which re- 
duces Eq. ([TJ to the Jaynes-Cummings Hamiltonian. As a con- 
sequence, the ground state of Eq. ([T]) is not a simple vacuum 
anymore as in the Jaynes-Cumming model. Instead, it con- 
tains finite average photon numbers, and shows non-classical 
properties such as squeezing and entanglement ITSl [T6l . To 
see this, let us examine the ground-state wavefunction more 
closely: Approximate expressions for the nearly-degenerate 
ground states have been derived in Ref. 17 (see also Ref.[T6]l. 
Here we take a different approach and explore the parity sym- 
metry in the Rabi Hamiltonian, which is important to under- 
stand the effect of photon hopping. The Hamiltonian ([T]l com- 
mutes with the "parity" operator 11; = exp(— i7ra|^ai)CT,f , and 
thus the Hilbert space is classified into subspaces £^ of ± 
parity. Within each subspace Sf, the Hamiltonian can be 
described in effect by a single bosonic operator, bi = fljcrf : 
^cQED ^ ^± ^ffO^ffl ^ith Hf = UJoibl - A/wo)(6. - 
A/wo) - AVwo and = f cos(7r6j6i) Oil. H° is sim- 
ply a displaced harmonic oscillator and the ground state is a 
coherent state jA/wo)^. For X/ojq ^ 1 (regardless of il). 
Hi can be treated perturbatively and shifts the energies of 
|A/ajo)b^ relatively by an exponentially small amount A = 
^e~^(-^/"'')^ Now, back in the {a^, erf }-basis, the nearly de- 
generate ground states \X/ujq)^ are expressed as 

|0),^^(|A/c.o)J+),-|-A/c.o),|-),), (3a) 

|l).^^(|AM)J+). + |-A/c.o),|->,), (3b) 

where \a)^ (a e C) is the eigenstate (coherent state) of 
and |±)^ are the eigenstates of af. In short, these two ground 
states, |0),- and residing in distinct parity subspaces are 
nearly degenerate with an energy splitting of 2A, separated 
far from higher-energy states by an energy gap ujq. 

Effective model: a transverse-field Ising chain — Let us 
now investigate the whole chain described by the Hamiltonian 
Typically J ^ ujq, and we are mainly interested in the 
low-lying excitations, well below uq. In this limit, each cQED 
remains within the subspace spanned by the states |0)j and 
in Eq. (j3|l and can be regarded as a pseudo-spin: 
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where rf — |0)j (0| — |1) ^ (1| and the energy splitting A plays 
the role of Zeeman field. Hopping of a photon into or out 
of a cavity changes the parity of its state, or more explicitly 
Qi |0)j = A/wo and a; |1) ■ = X/ojq |0)j. Based on these 



observation, we can identify and a] as a pseudo-spin-flip 
operator rf and aial_^_^ as Ising interaction t[t[_^_^. That is, 
the photon-hopping part of the Hamiltonian becomes 
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with Jeff = 2J{X/ll!o)'^ . The effective Ising interaction 
strength, J^s, is renormalized with respect to J by the fac- 
tor {X/luq)^ because the field part of the pseudo-spin states in 
Eq. (j3|l is a coherent state with amplitudes X/ujq and the field- 
field interaction between resonators is proportional to the am- 
plitudes of the resonator fields. 

Putting both terms in Eqs. Q and (jSj together, the low- 
energy effective Hamiltonian for the cQED chain becomes the 
so-called transverse-field Ising model (TFIM), 
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The TFIM exhibits a quantum phase transition between the 
magnetically ordered phase for A < J^s and the quantum 
paramagnet phase for A > J^s ESI- The former is partic- 
ularly interesting for our purposes. For A = 0, iJising has 
two degenerate ground states, |=>) = Hi and |<^=) = 
rij where |— and |<— ) • are eigenstates of rf . For 

A > (yet A < Jeff), rf tends to flip the pseudo-spins, 
^ causes tunneling between |=^) and \-^) via 

soliton propagation, and hence the true eigenstates become 
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However, as the tunneling involves N spins, the tunneling am- 
plitude is exponentially suppressed with the system size N. In 
other words, |5'o) and are nearly degenerate with energy 
splitting, S ^ exp(— A^/^) with ^ being the correlation length 
of the Ising chain, exponentially small in system size N. Both 
are separated from the continuum of excitations by the energy 
gap Jofi. 

The two states Iv&o) and in Eq. (j7|i have non-local 
combinations of many pseudo-spins and are widely known as 
Greenberger-Horne-Zeilinger (GHZ) states f25\. Moreover, 
by expressing them in the original {a;, erf } -basis 
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with s = or 1, one can see that they involve high degree 
of non-local entanglement between cavity fields and qubits. 
They are thus Schrddinger cat states at a truly large scale 
while many theoretically proposed or experimentally demon- 
strated Schrodinger cat states |26 , 27 1 contain merely a single 
radiation field. Below we illustrate that the two states in ([8]) 
are protected against local fluctuations and constitute a non- 
local qubit Ii20il . 
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Effective model: a Majorana chain — ID TFIM discussed 
above is equivalent to a chain of Majorana fermions | 28ll29l . 
The latter has attracted great interest because it permits local- 
ized Majorana modes that can be used for topologically pro- 
tected quantum computation Il28lj30ll . A very recent exper- 
iment 131] suggests that the Majorana chain can be realized 
in a solid-state system, and intensive efforts are made in this 
direction 1.32 J . 

Here we re-express the two nearly degenerate states in 
Eq. (j7]i or ([8]l in terms of localized Majorana fermions, and 
later discuss an experimentally feasible way of probing such 
Majorana fermions. The equivalence between the TFIM and 
the Majorana chain can be seen through a Jordan- Wigner 



transformation f33l: cj = t. 



-rf) with T, 



irf). The operators Ci and describe Dirac fermions and sat- 



isfy {ci,cj} = 5ij and {ci,Cj} = 0. The Dirac fermion op- 
erators are further represented with self-conjugate Majorana 



operators, 72^-1 = c, 
(|6]l is then reduced to 



and 72-^ ^ i{c\ - c,). The TFIM 



Majorana 



N 



N-l 



A 72j-l72i + -'"eft l2il2i+l 



(9) 



At A = 0, the Majoranas at the two ends, 71 and 72Ar, in the 
chain does not appear in the Hamiltonian, which implies the 
existence of two degenerate ground states. These are nothing 
but and |<;=) in Eq. (j7|i. For finite A, the two states |^) 
and |<^=:) are mixed linearly into |^'o) and j^I'i) in Eq. (|7j due 
to the tunneling between two Majorana modes 71 and 72Ar, 
and the degeneracy is lifted. Since the tunneling is through 
the whole chain, the energy splitting 5 is exponentially small 
(as long as A < Jeff). One can check that (71 +i72Ar) l^'o) = 
and (71 + i72w) = 2 |*o), which means that Iv^i) 
has one more fermion than I^Po) or equivalently that |^'o) and 
l^*!) have different fermion parities. 

Here we emphasize that the two Majoranas localized at the 
ends of the Majorana chain are actually non-local in the phys- 
ical chain, i.e., the cQED chain or the Ising chain f2T|: The 
Majorana operators is represented in terms of and as 
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and 72iv involves the string operator OjLiC^'''/)- This im- 
plies that the two nearly-degenerate ground states 1^*0) and 
l^&i) are not protected topologically against local noise even 
though mathematically they correspond to two distinct Majo- 
rana modes. It is in stark contrast to the case where the two 
Majorana modes at the ends of ap-wave superconducting wire 
are topologically protected. However, we will see below that 
the two states I^Po) and are vulnerable only to a certain 
type of local noise and there exist realistic systems with such 
type of local noise significantly suppressed. 

Noise — It is evident from the expression in Eq. (|8]l that 
the non-local spin qubits are prone to the local noise in erf of 



the physical qubits and the one in a; + a\ of the resonators. 
The states are intrinsically robust against the a\ and af noise 
since ^(vE'ilcr?''^ I^o)^ ^ e^i^/^a)'^ ^ which is reminiscent of 
the Franck-Condon effect. The + a\ noise affects only the 
resonator at the end of the chain (which is usually connected 
external microwave environment for measurement), and can 
be easily avoided by replacing it by a high-Q resonator The 
problem with erf noise can be circumvented, for example, by 
using Fluxonium for qubits. Fluxonium is known to have 
anisotropic noise characteristics with erf being the dominant 
noises and the erf noise ignorable fTSl. 

What about the effect of inhomogeneity in system parame- 
ters? Above we have assumed ^ ^nd J of each circuit- 
QED to be homogeneous. Deviations in u^, 17, A lead to fluc- 
tuations in A. The inter-cavity coupling strength, J, can also 
be varied from cavity to cavity, which leads to inhomogeneous 
TFIM, 
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This Hamiltonian still conserves the parity symmetry, V = 
riiLi '''1 which is respected by the degenerate ground states. 
Therefore, the ground states will be robust to small fluctua- 
tions in Aj and J, . We thus conclude that the nearly degener- 
ate ground states j^^o) and can be kept well protected by 
a careful design of the physical qubits in the system. 

Detection and control — In this section, we suggest a 
scheme to control and measure the non-local spin qubit. It 
can be also interpreted as detecting the Majorana bound states. 
Our proposal consists only of an additional empty resonator 
coupled to the resonator at the end of the circuit-QED chain. 
Consider a resonator with a frequency, oj^, capacitively cou- 
pled to A^th cavity, so that we have 



where ajv represents the field operator of A^th cavity, 
the field operator of the detection cavity. As shown ear 
A^th cavity's creation and annihilation operators are 
lent to X/ujqt^ for the A^th effective spin. Moreover, 
non-local spin qubits, rf is equivalent to S*^ = |4'o) 
|*i> (*o| for any i as rf l^-,) = l^-i^,) (s = 0, 1). 
fore, assuming that Jd = JdA/wo <C Jofi, the low- 
effective Hamiltonian (|6]) combined with the detection 
tonian ( 12 1 leads again to the Rabi Hamiltonian 



(12) 
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iJRabi =2^' + JdS^'iad + 4) + LUda^ad (13) 

Here we can make the rotating wave approximation, then the 
Hamiltonian reduces to the Jaynes-Cummings Hamiltonian. 
Therefore, by just adding an empty resonator at one end of 
the circuit-QED array, we can realize a circuit-QED Hamilto- 
nian for the non-local spin qubit. It allows us to tap into the 
standard techniques available for the circuit-QED to control 
and measure the non-local spin qubit. For example, since the 
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FIG. 2. (a) Energy diagram for the circuit-QED Hamiltonian l[TJ as 
a function of A/ctJo- (b) Plot of ^i{l\ai |0)^. We can conclude that 
A > 2ajo is required for our model to be valid because the transverse 
field A almost vanishes and the identification of photon annihilation 
operator as a spin flip operator, ^ = rf , is justified. 



detuning between the detection cavity frequency and the non- 
local spin qubit splitting, = Ud ~ S is large compared to 
Jd, it is in the dispersive regime where the cavity frequency 
pulling by the non-local spin qubit is 6ujd ~ 2J^/Ad ll23l . 
This can be experimentally measured since one can have 
Jd ~ lO^^oJa which achieves the standard strong-coupling 
regime for the circuit QED fS). 

Experimental feasibility — Finally we examine the exper- 
imental feasibility of the ideas explained above, estimating 
possible values of physical parameters of the system. Two 
requirements must be satisfied: First, the two ground states of 
each cQED in the system must be nearly degenerate and well 
separated from higher excitations. In Fig. Q (a) are plotted 
the energies of individual circuit-QED Hamiltonian ([T]l in the 
resonant case (ujq — Q). Figure |2](b) plots ^i{l\ai |0)j to 
illustrate how good (its value close to 1) the approximation 
Ui = X/ujqt[ is. One can see that A ^ 2ujq suffices for the 
requirement. Second, the system should be in the magneti- 
cally ordered phase (in terms of the effective TFIM), A < Jeff 
or equivalently fiexp [— 2(A/ajo)^] < 4J(A/wo)^- This re- 
quirement is satisfied provided that J > IO^^uq. The de- 
sired coupling strength, A > 2u!o, seems achievable for the 
Fluxonium coupled inductively to the superconducting res- 
onator |18|. Moreover, J > 10~^a;o is also realistic for the 



superconducting resonators, with J in the range of a few MHz. 

Conclusion - We have found several intriguing properties 
of the two nearly degenerate ground states of a chain of cou- 
pled circuit-QED systems in the ultrastrong coupling regime. 
The ground states are Schrodinger cat states at a truly large 
scale, and are mathematically equivalent to Majorana bound 
states. With a suitable design of the system, they are protected 
against local fluctuations, and may be probed and manipulated 
coherently by attaching an extra empty resonator. 

Finishing this work, we have noticed a closely related 
preprint ll34l . While they focus on the phase transition of the 
circuit-QED chain, we are mainly concerned about the quan- 
tum properties of the nearly degenerated ground states on one 
side of the phase transition. In this respect, both works are 
complementary to each other. 
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